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ABSTRACT 

We construct analytically stationary global configurations for both aligned and log- 
arithmic spiral coplanar magnetohydrodynamic (MHD) perturbations in an axisym- 
mctric background MHD disc with a power-law surface mass density Eo cc r~ a , a 
coplanar azimuthal magnetic field Bo oc r -7 , a consistent self-gravity and a power- 
law rotation curve vq oc r~@ where vq is the linear azimuthal gas rotation speed. The 
barotropic equation of state II oc E" is adopted for both MHD background equilibrium 
and coplanar MHD perturbations where n is the vertically integrated pressure and n 
is the barotropic index. For a scale-free background MHD equilibrium, a relation ex- 
ists among a, (3, 7 and n such that only one parameter (e.g., 0) is independent. For a 
linear axisymmetric stability analysis, we provide global criteria in various parameter 
regimes. For nonaxisymmetric aligned and logarithmic spiral cases, two branches of 
perturbation modes (i.e., fast and slow MHD density waves) can be derived once j3 is 
specified. To complement the magnetized singular isothermal disc (MSID) analysis of 
Lou, we extend the analysis to a wider range of — 1/4 < (3 < 1/2. As an example of 
illustration, we discuss specifically the (3 = 1/4 case when the background magnetic 
field is force-free. Angular momentum conservation for coplanar MHD perturbations 
and other relevant aspects of our approach are discussed. 

Key words: MHD — ISM: magnetic fields — stars: formation — galaxies: kinematics 
and dynamics — galaxies: spiral — galaxies: structure. 



1 INTRODUCTION 

Magnetic fields are ubiquitous in diverse astrophysical settings on various scales ranging from proto-stellar discs, young 
stellar objects (YSOs), microquasars, quasars, galaxies and grandiose astrophysical jets to clusters of galaxies. In some cases, 
magnetic fields are relatively weak so that dynamical processes are almost unaffected by their presence. However, there do 
exist numerous cases where magnetic fields are necessary and important for both dynamics and diagnostics, especially in 
spiral galaxies and accretion disc systems (e.g., Sofue et al. 1986; Beck et al. 1996; Balbus & Hawley 1998; Tagger & Pellat 
1999; Widrow 2002). In general, it is challenging to model magnetic fields realistically because of their complexities (Sofue 
et al. 1986; Kronberg 1994; Caunt & Tagger 2001; Widrow 2002). In the problems involving protostar formation and disc 
galaxies, Shu & Li (1997) presumed the so-called 'isopedic' magnetic field configuration where the mass-to-magnetic flux 
ratio in a razor thin disc remains constant. Using a singular isothermal disc (SID) model (Mestel 1963) that is isopedically 
magnetized, Shu et al. (2000) studied stationary coplanar perturbations and proposed that these global modes are bifurcations 
to either secularly or dynamically unstable configurations. Along a separate yet complementary line, Lou (2002) studied an 
azimuthally magnetized singular isothermal disc (MSID) and derived two different global stationary magnetohydrodynamic 
(MHD) perturbation modes. Based on the MSID model, Lou (2002) explored the manifestation of interlaced optical and 
magnetic field spiral arms in the outer portion of a disc with a nearly flat rotation curve such as the case of the nearby spiral 
galaxy NGC 6946 (Beck & Hoernes 1996; Fan & Lou 1996; Lou & Fan 1998a, 2002; Frick et al. 2000, 2001; Lou et al. 2002). 

In contexts of spiral galaxies, it is natural and important to consider a composite disc system consisting of gravitationally 
coupled gaseous and stellar discs. This is because various physical processes on different scales occur in the gaseous disc yet 
large-scale gas dynamics and environment are significantly affected by large-scale structures in the stellar disc (Lou & Fan 
1998b). To construct global coplanar perturbation structures in a systematic manner, we started with a composite SID system 
of two fluid discs without magnetic field (Lou & Shen 2003) to generalize the work of Shu et al. (2000) for a single SID. In 
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terms of the global axisymmetric stability for such a composite SID system, we provided a more straightforward D— criterion 
(Shen & Lou 2003) in contrast to the axisymmetric stability criteria of Elmegreen (1995) and Jog (1996). To be more general 
than a background SID profile, we recently constructed global stationary perturbation patterns in a composite system of two 
fluid scale-free discs without magnetic field (Shen & Lou 2004a, b) to further generalize the work of Syer & Tremaine (1996). 
With proper adaptations, these global stationary perturbation solutions may be utilized to model large-scale structures spiral 
galaxies. 

While a scale-free disc has no characteristic scales, one may impose boundary conditions to describe a modified scale-free 
disc, e.g., by cutting a central hole in a disc (e.g., Zang 1976; Evans & Read 1998a, b). With such boundary conditions, discrete 
eigenfunctions for perturbations may be constructed as growing normal modes. On the other hand, by specifying a phase 
relation for a postulated reflection of spiral waves from the origin r = 0, Goodman & Evans (1999) could also define discrete 
normal modes even for an unmodified gaseous SID. Shu et al. (2000) speculated that the swing amplification process (e.g., 
Goldreich & Lynden-Bell 1965; Julian & Toomre 1966; Toomre 1977; Fan & Lou 1997) across corotation allows a continuum 
of normal modes (Lynden-Bell & Lemos 1993) while proper 'boundary conditions' may select from this continuum a discrete 
spectrum of unstable normal modes; they also suggested that zero- frequency (stationary) disturbances do signal the onset of 
instability. 

On the basis of the work of Shu et al. (2000) for a single isopedically magnetized SID and of Lou & Shen (2003) for 
a composite SID system without magnetic field, Lou & Wu (2004) analyzed coplanar MHD perturbation structures in a 
composite SID system with one of the disc being isopedically magnetized (Shu & Li 1997). As X-ray emitting hot gases can 
freely flow along open magnetic field lines anchored at the gaseous disc, global stationary solutions constructed by Lou & 
Wu (2004) set an important preliminary stage for modelling spiral galactic MHD winds. In particular, if activities of star 
formation, starbursts and supernovae etc. are mainly responsible for producing X-ray emitting hot gases streaming out of the 
galactic plane from both sides, we would expect stronger winds from the circumnuclear starburst region (e.g. Lou et al. 2001) 
and along spiral arms than those winds from the remainder of the galactic disc plane. 

Both SID and MSID models offer valuable physical insight for disc dynamics and belong to a wider family of scale-free 
discs (Syer & Tremaine 1996) where physical variables in the axisymmetric background equilibrium scale as powers of radius 
r (e.g., the rotation curve oc r~ 13 ). It turns out that coplanar perturbations in such scale-free discs in a proper range of 
(5 can be treated globally and analytically (Lemos et al. 1991; Syer & Tremaine 1996) without invoking the usual WKBJ 
or tight-winding approximation to solve the Poisson equation for density wave perturbations* (Lin & Shu 1964; Binney & 
Tremaine 1987; Bertin & Lin 1996). In other words, it becomes feasible to include azimuthal magnetic fields in a scale-free 
disc and to construct global coplanar MHD perturbation configurations that are stationary in an inertial frame of reference. 
Such analytical solutions, still idealized and limited, are making important further steps and are extremely valuable for 
bench-marking numerical codes and for initializing numerical MHD simulations. 

For an azimuthally magnetized scale-free disc with an infinitesimal thickness, the axisymmetric background equilibrium 
without radial flow is characterized by following features: a surface mass density Eo oc r~ a , a rotation curve vo oc r~^, an 
azimuthal magnetic field Bo oc r~ 7 and a barotropic equation of state n = KT, n where n is the vertically integrated gas 
pressure with constant K > for a warm disc and n > as the barotropic index. To maintain a radial force balance at 
all radii, there exists an explicit relationship among a, f3, 7 and n such that only one parameter is independent. As will be 
seen in Section 2, it is convenient to use (3 as this independent parameter (Syer & Tremaine 1996; Shen & Lou 2004b). In 
order to satisfy relevant constraints, the prescribed range of f3 turns out to be —1/4 < (3 < 1/2 that includes the previous 
special case of /? = (i.e., the MSID case) studied by Lou (2002), Lou & Zou (2004a) and Lou & Wu (2004). In Section 2, we 
introduce coplanar MHD perturbations and present the linearized perturbation equations. We construct and analyze global 
stationary aligned and logarithmic spiral perturbation configurations in Section 3 and as an example of illustration, we discuss 
specifically the case of (3 = 1/4 with the background azimuthal magnetic field being force-free. In Section 4, we derive phase 
relationships between perturbation enhancements of surface mass density and magnetic field for both aligned and logarithmic 
spiral cases, and discuss the problem of angular momentum conservation. Finally, we summarize our results in Section 5. For 
the convenience of reference, relevant technical details are collected in the Appendices. 



* By numerically solving the perturbed Poisson equation using the Fourier-Bessel transform, Riidiger & Kitchatinov (2000) approached 
the linear stability of the disc in terms of an eigenvalue problem and could also treat the density wave problem globally. Their formulation 
differs from ours in two major aspects. Firstly, they include a central mass that mainly determines the rotation curve (i.e., an approximate 
Keplerian disc rotation for a disk mass being much less than the central mass). In our model, the disc rotation curve is completely 
controlled by the disc self-gravity and our approach is semi-analytical. Secondly, for axisymmetric instabilities, they only revealed ring 
fragmentation instabilities when the rotational Mach number becomes sufficiently high (which is similar to our result). However, our 
model analysis further indicates another type of collapse instabilities (i.e., magneto-rotational Jeans instability) in the perturbation 
regime of large radial spatial scales (see Section 3.2.1 for more details). 
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2 FORMULATION OF THE PROBLEM 

In the magnetofluid approximation, the disc is taken to be razor-thin (i.e., we use vertically integrated magnetofluid equations 
and neglect vertical derivatives of physical variables) and large-scale stationary aligned and spiral coplanar disturbances 
develop in a background MHD rotational equilibrium of axisymmetry (Syer & Tremaine 1996; Lou 2002; Lou & Zou 2004a; 
Lou & Wu 2004). The background magnetic field is taken to be azimuthal to avoid the magnetic field winding dilemma (Lou 
& Fan 1998a). For the sake of simplicity at this stage, non-ideal effects such as viscosity, resistivity and thermal diffusion etc. 
are ignored for large-scale perturbations. 



2.1 Basic Coplanar MHD Equations in a Cylindrical Disc Geometry 

Using cylindrical coordinates (r,6,z), we have the following two-dimensional nonlinear MHD equations for a razor-thin disc 
geometry: the mass conservation equation 

as i djrZu) i a(sj) _ 

dt r dr + r 2 89 ' { ' 

where E is the surface mass density, u is the radial bulk flow velocity, j = rv is the specific angular momentum in the vertical 
z-direction and v is the azimuthal linear velocity; the radial component of the momentum equation 
du +u du + j_du f_ _ _ l_an d(j> T 1 fdzBg\d(rB g ) dB r 



dt dr r 2 88 r 3 E dr dr E / 47rr 



(2) 



dr d9 

where n is the vertically integrated gas pressure (sometimes referred to as two-dimensional pressure), 4>t is the total grav- 
itational potential including that of an axisymmetric dark matter halo in contexts of disc galaxies, and Be and B r are the 
azimuthal and radial components of the coplanar magnetic field B, respectively; the azimuthal component of the momentum 
equation 

~d(rB g ) dB r ' 



dr d6 



(3) 



dt dr r 2 d9 E d6 39 E J 4tt 
the Poisson integral equation 

F** = -af**l [r > 2 + r 2 -2r>rUi,-e)\V 2 ' (4) 
where F = 4>/4>t is defined as the ratio of the potential arising from the disc to that arising from the entire system including 
a massive dark matter halo that is presumed not to respond to coplanar perturbations in the disc plane (Syer & Tremaine 
1996; Shu & Li 1997; Shu et al. 2000; Lou & Shen 2003; Shen & Lou 2004a, b; Lou & Zou 2004a) with F = 1 for a full disc 
and < F < 1 for partial discs; the divergence- free condition of magnetic field B = (B r , Bg, 0) 

+ 

the radial component of the magnetic induction equation 

dB r Id,. . . 

-W = rde {uBe ~ vBr) ' (6) 

the azimuthal component of the magnetic field induction equation 

(7) 

and the barotropic equation of state 

n = KTT , (8) 

where constants K > and n > 0. Among the three equations JBJ and Q, we can freely choose two independent ones. 
By barotropic equation of state Q, the sound speed a is defined by 

2 _ dJlo n _i , . 

a = — = nKT, , (9) 

CLIjO 

where subscript indicates the background equilibrium. An isothermal sound speed corresponds to a barotropic index n = 1. 



2.2 Properties of an Axisymmetric Equilibrium 

We now proceed to derive properties of an axisymmetric rotational MHD equilibrium characterized by physical variables 
associated with a subscript 0. In our notations, such a background equilibrium has the following power-law scalings: surface 
mass density of Eo (X r~ a , a radial velocity profile uo = and jo = rvo oc r 1- ' 3 , a purely azimuthal magnetic field with Br = 
and Bg = Bo oc r Substitutions of these power-law scalings into equations — (0 lead to the following condition for the 
radial force balance, namely 

^+aa 2 -(l- 7 )C7i = rM , (10) 
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where the Alfven speed Ca in a thin disc is defined by 

c ^ikI B2 ° dz 

and the Poisson integral gives 

Fcj>T = -2irGrZ P {a) , (12) 
with the numerical factor Vo(ot) explicitly defined by 
-n r y_ r(-q/2 + l)r(q/2-l/2) 

V ° {a) = 2 r(-a/2 + 3/2)r(a/2) (13) 
and r(...) being the standard gamma function. This expression 1131 can also be included in a more general form of Vm(f3) as 
in expression 1421 derived later with the limiting result of 2/3'Po(f3) — » 1 as (3 — > 0. 

To satisfy the radial force balance equation (I1UI at all radii (i.e., scale-free) would necessarily require^ 

2P = a(n - 1) = 27 - a = a - 1 , (14) 
which immediately leads to the explicit expressions of indices a, 7 and n all in terms of [3 

In this formulation, <^ remains finite only if 1 < a < 2 which implies < (3 < 1/2. For the total gravity force arising 
from the equilibrium surface mass densities to be finite, a larger /3— range [3 6 ( — 1/2, 1/2) is allowed. The (3— range is also 
constrained by —1/4 < j3 < 1/2 where the first inequality is required for n > and the second inequality is imposed such 
that the central point mass will not diverge. By these considerations, the plausible range of (3 falls in ( — 1/4, 1/2) (see Syer & 
Tremaine 1996 and Shen & Lou 2004a, b). 

According to equations (11 Ot and 1150 . we can explicitly write all physical variables of the background equilibrium state 
in power-law scalings of r in terms of parameter f3, namely 

FA 2 [1 + D 2 -(l-4/3)g 2 /(2 + 4/3)] 2 A 2 2 _ 2 2 q 2 A 2 

^0 - n _„, na ~. r-TZM 1 a - /-, , na \..oei > C A = q a - 



2-KG(2f3Vo)r 1+20 ' (1 + 2/3)^' A ~ H (1 + 2/3)^' 

«o = ^, j =rv = ^, n = ^ = ^, K =[(2fi/r)d(r 2 n)/dr] 1/2 = [2(l-/3)] 1/2 n. 

By introducing a reference radius 7?, the constant A is actually related to a scaled sound speed a(l + 20) 1 ^ 2 R fi , the constant 
D = i>o/[a(l + 2/3) 1 ' /2 ] is a scaled rotational Mach number, the constant q is the ratio of the Alfven speed Ca to the sound speed 
a, Q and k are the equilibrium angular speed and the epicyclic frequency, respectively (Lou 2002; Shen & Lou 2004a, b). The 
first expression of equation 11161 puts an upper limit on the magnetic field strength, namely, 1 + D 2 — (1 — 4f3)q 2 /(2 + 4/3) > 0. 
Alternatively, once a magnetic field parameter q 2 is chosen, the scaled rotational Mach number D must meet the following 
physical requirement 

">%^- 

for a positive equilibrium surface mass density Eo. Note that for 1/4 < f3 < 1/2, inequality 1171 is automatically satisfied 
without actually restricting q 2 and the relevant physical requirement is simply D 2 > 0. Not all stationary solutions of D 2 can 
satisfy condition 1171 . This then implies limits on allowed magnetic field strength for a known background rotational MHD 
equilibrium. 

It is also important to note that for the axisymmetric background rotational MHD equilibrium under consideration, 
the magnetic force arising from the azimuthal magnetic field is radially inward when —1/4 < (3 < 1/4, radially outward 
when 1/4 < f3 < 1/2 and is zero when (3 — 1/4 [i.e., the azimuthal magnetic field is force-free (Lou & Fan 1998a) and the 
equilibrium density distribution reduces to that of the single unmagnetized disc case (Lemos et al. 1991; Syer & Tremaine 
1996)]. These different possibilities can be physically understood in terms of the competition between the magnetic pressure 
and tension forces. In order to see this more specifically for the background, we can write -(47rEr)~ 1 / dzB Q d{rB ) / dr = 
— E -1 [(<9/cV) f Bq(1z/(Stv) + f B 2 dz/(47vr)] referring to equation The first term in the square brackets stands for the 
magnetic pressure force and the second term stands for the magnetic tension force for a purely azimuthal magnetic field 
Bo(r) in the cylindrical geometry. As Bo scales r -7 , we have (d/dr) f Bgdz/^n) + f B 2 dz/(4Trr) = (—7+ 1) f B 2 dz/(4nr). 
Therefore for < 7 < 1 or equivalently —1/4 < [3 < 1/4, the magnetic tension force dominates and the total magnetic force is 
radially inward; for 1 < 7 < 3/2 or equivalently 1/4 < [3 < 1/2, the magnetic pressure force dominates and the total magnetic 
force is radially outward. The two magnetic forces balance each other for 7 = 1 (or (3 = 1/4). 



t The magnetic field must obey the scale- free requirement H14I in general and for the special case of (3 = 1/4, we have 7 = 1 for the 
equilibrium magnetic field to be force-free (e.g., Low &: Lou 1990). The non- force-free MSID case studied by Lou (2002) has 7 = 1/2, 
f} = 0, a = 1 and n = 1. 
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2.3 Coplanar MHD Perturbations in the Disc Plane 

On the basis of the full MHD equations Q — J7J), we readily derive the linearized coplanar MHD perturbation equations as 
9Ei 1 9(rE ui 



dt 

du\ 
~~dt 

dji 



+ n 



dr 

du\ 2Qji 
~~d6 r~ 



de 



o 



d_ 

dr 
d 



a 2 S! 
So 



+ < 



(l-4/3)CiE! 



dt + 2n 1+ de ~ de 



a 2 Ei 



= -G <l> dip 



So 

Si(r',^,t) r'dr' 



S 



2S r X 
8(rBo) 



dzBo 



47rr 



d(rbg) db r 



So / 



dzb e d(rBo) 
4irr dr 



4tt 



dr 



[r' 2 + r 2 - 2r'r cos(^ - #)] 1/2 



d(rb r 
dr 

db r 
~dt = 



db e 



+ 
a 

ld_ 

r d9 

t) 







(tti-Bo — rQ,b r 
(uiBq — rQ.br 



(18) 
(19) 
(20) 
(21) 
(22) 
(23) 
(24) 



dt dr 

where we use subscript 1 to denote associations with small disturbances in physical variables and b = (b r ,bg, 0) stands for the 
coplanar magnetic field perturbation. As the background rotational MHD equilibrium is stationary and axisymmetric, these 
perturbed physical variables can be decomposed in terms of Fourier harmonics with the periodic dependence exp(iiut — im,e) 
where ui is the angular frequency and m is an integer for azimuthal variations. More specifically, we write 

Si = S(r) exp(kji — im#) , u\ — U(r) exp(kjt — im(9), ji — J(r) exp(kji — imff) , 

4>i = V(r) exp(kjt — \m.e) , b r — R(r) exp(kjt — im(9) , be = Z(r) exp(kji — imff) , 

where S(r), U(r), J(r), V(r), R(r) and Z(r) are all radial variations of the corresponding perturbed physical variables and 
can be complex in general for possible radial oscillations. Without loss of generality, we take m > in our analysis. 

With Fourier decomposition 12511 . it is then straightforward to cast coplanar MHD perturbation equations 1181 1 — 1241 1 
into the following forms of 

imSo 



(25) 



Ulj - mQ)S + --^-(rEoU) 
r dr 



-J = 



i(u> — mQ)U 



2Q.J 
r 



dr 



(1 - 40)C\S 
2S r 



1 



dzZ d(rB ) 
Airr dr 



1 



dzBo 
4nr 



d(rZ) 



dr 



+ imR 



where $ = a 2 5/So + V is a short-hand notation, and 



rK I 
i(w - mQ)J + -—-U = 1771$ + — 

2l I 2-iQ 



V{r) =-G (p dip 
d(rR) 



f 

Jo 



dzR d{rB ) 
4ir dr 

S(r') cos(rmp)r'dr' 
(r' 2 + r 2 — 2rr' cos xp) 1 / 2 ' 



dr 



- imZ = , 



i(cj-mQ)R+^^-U = 



™ z = iirSIR) - £-(B U) 



(26) 
(27) 

(28) 
(29) 
(30) 
(31) 
(32) 



dr dr 

where equation IH2t can be derived by combining equations l-'ilH and l|31lL 

We now rearrange the time-dependent coplanar MHD perturbation equations by taking m > 1 from equations l|26|l — l|32|l 
and the special case of m = will be analyzed in details at the end of this subsection. 
From equations 1301 and 131H . we readily obtain 
i d(rR) „ mB U 



m dr 



R : 



r{u> — mQ.) 



(33) 
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A substitution of equation 1331 into the radial and azimuthal components of the momentum equation I27H and 128H yields 



i((jj - 
and 



mQ)U 



2Q.J 



i(u - mfl)J + —U 



dr 



im<E> 



(1 - 4/3)CjS 
2S r 



mCi(l - 4/3) 
2r 



U 



dr- 



+ 



12/3 d 2/3(1 + 4/3) 

"T" o 



2r dr 



U 



(34) 



(35) 



respectively. 

In order to construct global solutions without the WKBJ approximation, we are mainly interested in stationary config- 
urations with zero pattern speed u>/m = 0, also referred to as neutral ui = modes (e.g., Syer & Tremaine 1996; Shu et al. 
2000; Lou 2002; Shen & Lou 2004a, b; Lou & Zou 2004a). By setting to = and taking m > 1, equations J33J|, fHJ and 

13511 can be readily reduced to 



„^ 1 d , _ . Trx mSn i 
mQ,S + - — (rE iU) + — 5— J ; 
r or r l 







(36) 



^_ 9 (B iU\ 
dr \ mft J 



R 



B U 



(37) 



n . rf 29. J 
mQ.iU H = — 



9$ 

9r 



(1 



-4/3)CiS 



r 2 



2E r 



+ 



12/3 9 2/3(1+4/3) 



2r 9r 



if/ 



(38) 



mOZ+^it/ 



Ci(l-4/3) iC7 



(39) 



2r ' 

Equations 1361 through 139H together with equation 129H are the basic MHD perturbation equations for constructing 
stationary non-axisymmetric m > 1 configurations of both aligned and unaligned logarithmic spiral cases. Note that aligned 
and spiral stationary global solutions both involve propagations of fast and slow MHD density waves (Lou 2002). 

We next present the basic coplanar MHD perturbation equations for stationary axisymmetric m — configurations with 
or without radial propagations. Starting from equations I26L l|31|l . ll:-)2> . l|34|l and l|35|l by setting m = and u) 7^ 0, we obtain 



1 d 



r dr 
2ujQJ 



(rE i)7) 



9$ (1 
or 



4P)C\S 



2E r 



R = , 



Ll)Z ■ 



d_ 

dr 



(BoiU) 



dr- 



+ 



12/3 d 2/3(1 + 4/3) 



2r dr 



iU 



(40) 



This set of equations is applicable as we proceed to investigate stationary axisymmetric coplanar MHD perturbations (Lou 
& Zou 2004a, b). For coplanar hydrodynamic perturbations or coplanar MHD perturbations in an isopedically magnetized 
SID (Shu et al. 2000), there is no essential difference for changing the order of limiting procedure for m = and uj — > 0. 
However, for MSID with a coplanar azimuthal background magnetic field, the results would be different by changing the order 
of limiting procedure for m = and u — > (Lou 2002; Lou & Fan 2002; Lou & Zou 2004a). 



3 STATIONARY ALIGNED AND LOGARITHMIC SPIRAL CONFIGURATIONS 

To solve the Poisson equation connecting the perturbed surface mass density and the perturbed gravitational potential, we may 
consistently assume disturbances to carry proper scale-free forms (Lynden-Bell & Lemos 1993) for aligned and logarithmic 
spiral perturbations. For aligned perturbations, S(r) contains only an amplitude variation in r of perturbed surface mass 
density and does not involve phase variation in r so that the maximum density perturbations at different radii line up in 
the azimuth. For spiral perturbations in comparison, in addition to an amplitude variation in r, S(r) also involves a phase 
variation in r such that a spiral pattern emerges. 



3.1 Aligned Global Coplanar MHD Configurations 

For aligned global coplanar MHD perturbations, we select those perturbations that carry the same power-law variation in r 
of the background equilibrium, namely 

S = ar- 1 - 20 , V = ~2nGrST m (f3) , (41) 
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where a is a small constant coefficient and the numerical factor V m {P) is given explicitly by 

v _ r(m/2- / 3 + l/2)r(m/2 + /3) . . 

mKP) 2r(m/2-/3 + l)r(m/2+/? + l/2) ' 1 ' 

where -m/2 < 8 < (m + l)/2 (Qian 1992; Syer & Tremaine 1996; Shen & Lou 2004a, b). In the limit of (3 -» 0, we simply 
have P m — 1/m consistent with earlier results on SIDs and MSIDs (Shu et al. 2000; Lou 2002; Lou & Shen 2003; Lou & Zou 
2004a; Lou &i Wu 2004). In fact, a more general class of self-consistent potential-density pairs satisfying the Poisson integral 
1291 is also available; see footnote 3 in Shen & Lou (2004b). 



3.1.1 Aligned m = Case 

For aligned neutral modes of axisymmetry (m = 0), we can start from equation 140H by setting ui = 0; this implies U = R = 
and gives no constraints on Z. As S oc r _1_2/3 , the scale-free condition requires Z oc r -7 with 7 = 1/2 + 2/3 and J oc r 1- ' 3 . It 
turns out that such stationary coplanar MHD perturbations are simply alternative states for the axisymmetric background 
equilibrium with a proper rescaling factor (Shu et al. 2000; Lou 2002). As this is somewhat trivial, we now turn our attention 
mainly to cases of non-axisymmetric m > 1 stationary aligned coplanar MHD perturbations. 



3.1.2 Aligned m > 1 Cases 

For the potential-density pair S oc r _1_2/3 and V = —2-KCPmrS , we have U oc r - ' 3 and J oc r 1- ' 3 according to equations 
<5Et - O. It then follows that 

mfiS-^E i[/+^J = 0, 



, , /v. ma 2 / rEo „ A „ C 2 A (l — 48) , rT 

mSlJ + r(l - /?)Oi{7 = —f-2-KGV™ -1 rS+ AV of ^ — — it/ 

rLo V a / 2S2r 



(43) 



2 y milr fir 
Rearranging the first three equations of (0HJ, we obtain a set of three algebraic equations, namely 

a x S + 6ii(7 + ciJ = , e^S + b 2 i[7 + c 2 J = , a 3 5 + b 3 i(7 + c 3 J = , (44) 

where coefficients a%, &», Cj (i = 1, 2, 3) are explicitly defined by 

n , 3/3E mE 
ai = mil , 01 = , ci 



_ (l-4/3)Cj 
fl2 = 2S r ~ 2/3 



2nGV m (8) - 



[(1 - 40) 2 - 2m 2 \C 2 A _ 2D. 



&2 =™" + ' C2 = V («) 



03 



2 

27rGP m (/3) - — 
Eor 



2mQr 2 r 

b3= r (i-p)si- (i ~y )c ^ ■ c 3=™" 



For nontrivial solutions of {S 1 , i(7 , J}, the determinant of algebraic equations l|44|l must vanish, namely 

6i(o 2 c 3 - GE3C2) ci(a 2 &3 - CI3&2) n ,.m 

ai 71 r n 1 71 1 r~ = u • l 4t) i 

(o 2 c 3 - o 3 c 2 ) (& 2 c 3 - o 3 c 2 ) 

Using coefficient definitions 145H in condition 14611 . we obtain a more informative form of stationary dispersion relation 

a 2 - 2TrGP m T, r - Q 2 r 2 Q. 2 r 2 [68a 2 + (2 - 20)Q 2 r 2 - l28irGV m Z r - (1 - 4/3) Cj] 

(l/2-2/3)Ci-(l + 2/3)Sl 2 r 2 ~ [ra 2 Sl 2 r 2 - Ci(m 2 - 1/2 - 8/3 2 + 4/3)] [(1 + 2/3)fi 2 r 2 - Ci(l/2 - 20)} 

_ 40a 2 - 8/37rGP m S r + (1 - 4/3) Cj 
2[m 2 fi 2 r 2 - C\(m? - 1/2 - 8/3 2 + 4/3)] " 

Note that equation 147H takes a simpler form in the limit of — > (i.e., a flat rotation curve of a SID) and reduces to equation 
(3.1.9) of Lou (2002) as expected. 

By using background MHD equilibrium variables 1161 in condition 14711 . we obtain the stationary dispersion relation for 
global aligned m > 1 coplanar MHD configurations in the form of a quadratic equation in terms of y = D 2 , namely 

C 2 y 2 + C x y + Co = , (48) 
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where the three coefficients d (i = 1, 2, 3) are defined by 
" (3 - 4/3)m 2 - (1 - 4/3)(l - 6/3 + 4/3 2 ) , . 



Ci 



Co 



FCPm + m - (1-4/3) (2-/3) 



g -An(l-FCT> m ) , 



4/3 



4(1 + 2/3) 



2(1 + 2/3) 

2m 2 + (1 - 4/3) (1 - 2/3) 



(49) 



2m 2 + (1 - 4/3) (1 - 2/3) 
1 + 2/3 



FCP m + 1-4/3 



with auxiliary parameters explicitly defined by 



A m = m 2 + 4/3 - 4/3 2 , B m = (l + 2/3)(m 2 -2 + 2/3) , C = , W m = FCP m An + B m (50) 

2pFo 

which are exactly the same as those defined in Shen & Lou (2004b) for aligned cases of global perturbations in a purely 
hydrodynamic composite disc system. Here, the effect of coplanar magnetic field is associated with the q 2 parameter in 
definitions 1491 . The determinant of quadratic equation 1481 . A = C 2 — 4C2C0, depends on various parameters and can 
become negative under some circumstances (i.e., no physical solutions for y = D 2 ). For A > 0, the two real solutions of 
y = D 2 to equation l|48[l are 

-C\ + A 1/2 -Ci - A 1 / 2 

»i = — — and = ^c 2 — ■ (51) 

As a first check of necessary consistency, we consider the limiting case of vanishing magnetic field, that is, q — 0. Then, 
equation 1481 simply reduces to 

H m y 2 - A m (l - FCV m )y = , (52) 
which gives a non-trivial solution 

D 2 — y\= ^— FCVm)Am ^ (53) 

This D 2 solution 1531 is exactly the same result for the case of a single scale-free gas disc without magnetic field (see Syer & 
Tremaine 1996 and subsection 3.1 of Shen & Lou 2004b). 

As a second check of necessary consistency, we consider another limiting case of /3 — > which corresponds to the MSID 
case with a flat rotation curve studied by Lou (2002). One can readily show that equation 11481 reduces to equation (3.2.11) 
of Lou (2002) as expected. 

In contrast to the case of a single fluid disc studied by Syer & Tremaine (1996) without azimuthal magnetic fields, there 
are now two branches of D 2 solutions of algebraic equation 1481 in general. A D 2 solution is considered to be physical when it 
satisfies both the condition D 2 > and inequality l|17|l . When the magnetic field becomes sufficiently weak, among these two 
branches of D 2 solutions, y\ is the counterpart of the single disc case and 1/2 is additional due to the coupling between the 
surface mass density and magnetic field. We recall that in a composite system of one stellar disc and one gaseous disc coupled 
through the mutual gravitational interaction, there are also two branches of D 2 solutions (Lou & Fan 1998b, 2000; Lou & 
Shen 2003; Shen & Lou 2004b). In a composite stellar-gaseous disc system, the two different classes of perturbation modes 
correspond to either in-phase or out-of-phase of surface mass density perturbation enhancements (see also Lou & Wu 2004). 
In a single azimuthally magnetized disc, the two coplanar MHD modes yi and 1/2 will be distinguished by either in-phase or 
out-of-phase perturbation enhancements of the surface mass density and the azimuthal magnetic field in the WKBJ or tight- 
winding approximation. In the WKBJ regime, the y\ and yi branches correspond to stationary fast and slow MHD density 
waves, respectively (Fan & Lou 1996; Lou & Fan 1998a). Although the phase relationships for perturbation enhancements 
of the surface mass density and the azimuthal magnetic field become ±7r/2 for open spiral structures (Lou 2002; Lou & Fan 
2002), we shall still refer to y\ and 2/2 solutions as global stationary fast and slow MHD density waves, respectively. We shall 
come back to this discussion of 'phase relationship' in more details in Section 4. 



4 -r-vi-^u-^ fl FCV )a 2 

q 2(1 + 2/3) 1 ^l m )q 



1 + 2/3 



3.1.3 Force-Free Magnetized Discs with f3 — 1/4 

As an example of illustration, we now focus on the specific case of /3 = 1/4 with inequality l|17|l automatically satisfied. In 
such a magnetized disc system, the rotation curve scales as r _1//4 , the surface mass density scales as r _3//2 , the barotropic 
index is n — 4/3 and the azimuthal magnetic field scales as r _1 (i.e., the background equilibrium magnetic field is force-free). 
For this specific case, the determinant A of algebraic equation 1481 remains always positive for both full and partial discs (see 
proofs in Appendix A). As seen from D 2 solutions 1511 . the sign of C2 determines whether 1/1/1/2 is the upper/lower branch or 
reverse. For m > 2, we always have C2 > and therefore y\ and 1/2 remain to be the higher and lower branches, respectively. 
The situation becomes somewhat involved for the m = 1 case when the sign of C2 is dependent on the value of F parameter. 
We show in Fig. 1 the y = D 2 solution branches for m = 1 and m = 2 cases, respectively, versus the gravitational potential 
ratio F. 

For the m = 1 case, there exists a diverging point at F c — 3/[7CPi(l/4)] for yi = D 2 solution branch where C2 = Hi = 0. 
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(a) m = 1, /3 = 1/4, g = 0, 1 



(b) m = 2, /3 = 1/4, q = 0, 1 



Figure 1. Two solution branches for y = D 2 of algebraic equation 1481 for m = 1, 2, (3 = 1/4 and g = 0, 1. The j/i branch is plotted 
in solid lines and the j/2 branch in dash lines. Panel (a): For the m = 1 case, the relative locations of the two branches switch across 
the diverging point at F c = 0.6842. Physical solutions are those above the y = axis marked by the light dash line. For < F < F c , 
only the y2 branch (the dash lines) is physical. For physical solutions, the enhancement of magnetic field will raise both branches but 
with different rate. The y\ branch (the solid lines) is apparently influenced by the presence of the coplanar magnetic field, while the y2 
branch (the dash lines) is slightly shifted. Panel (b): For the m = 2 case, y\ and j/2 solutions remain always to be the upper and lower 
branches, respectively, and both branches are physical in the range of < F < 1. Both branches are raised by the enhancement of the 
coplanar magnetic field strength. The m > 2 cases are qualitatively similar to the m = 2 case shown in panel (b) here. 



The value of this critical F c is numerically determined to be F c ~ 0.6842. More specifically, we have C2 < for < F < F c , 
while we have C2 > for F > F c . This sign variation of C2 is indicated by the relative locations of 1/1 and 1/2 solution branches 
as shown in panel (a) of Fig. 1. In other words, for < F < F c , y 1 remains negative and the only physical solution is 2/2- On 
the other hand for F > F c , both y\ and 1/2 solution branches become physically plausible. 
For j3 = 1/4, the explicit form of algebraic equation 1481 becomes 

D 2 + |(1 - FCV m )m 2 q 2 = , (54) 

where auxiliary parameters P m , A m , B m , C and 7i m are all evaluated at (3 = 1/4. It is fairly straightforward to show that 
< CP m < 1 for all m > 1 when (5 = 1/4 (Shen & Lou 2004b). As < F < 1, we then have < 1 - FCV m < 1. 



FCP„ 



2 3 
m +- 



3 / 2 
2 {m 



D 4 



( 2FCP„ 
V 3 



1 )m 2 q 2 + 



Ik. 



FCPr, 



3.2 Logarithmic Spiral Configurations for Global Coplanar MHD Perturbations 

In our formulation, stationary surface mass density perturbations are characterized by S(r) exp(— im9). For aligned coplanar 
MHD perturbations, we took S oc r~ e where e is a real exponent. In Section 3.1, we have chosen e — a = 1 + 2/3 that carries 
the same power-law dependence as the equilibrium disc does. For a complex e, a perturbation pattern would appear spiral, 
that is, a logarithmic spiral in the form of S oc r" 34 ^' exp[— iSJ(e) lnr] where Sft(e) and Ss(e) are real and imaginary parts 
of e. To ensure the gravitational potential arising from this perturbed surface mass density being finite, we should require 
— m + 1 < K(e) < m + 2 (e.g., Qian 1992). In fact, there exists a more general class of self-consistent potential-density pairs 
as indicated in footnote 7 of Shen & Lou (2004b). However, for the analysis of axisymmetric stability problem in Section 
3.2.1, equation 1551 is used such that dispersion relation I57H derived later on is real on both sides (e.g., Lemos et al. 1991; 
Syer & Tremaine 1996; Shu et al. 2000). To be specific here, we consistently take the logarithmic spiral surface mass density 
perturbations and the resulting gravitational potentials as (Kalnajs 1971; Lemos et al. 1991; Shu et al. 2000; Lou 2002; Lou 
& Shen 2003; Shen & Lou 2004b; Lou & Zou 2004a) 

S =0T -3/2 exp(i£lnr) and V = -27rGrS7V m (£) , (55) 

where cr is a small constant, £ is a parameter related to the radial wavenumber and the Kalnajs function A/" m (£) is defined by 

_ 1 r(m/2 + ig/2 + l/4)r(m/2 - ig/2 + 1/4) 
m ^> 2r(m/2 + i£/2 + 3/4)r(m/2-i£/2 + 3/4) {ou > 

(Kalnajs 1971). As A/" m (£) is an even function of £, a consideration of £ > would suffice. In our convention of notations, 
(, > and £ < correspond to leading and trailing logarithmic spiral waves, respectively, and £ = kr relates the parameter 
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£ to the radial wavenumber k. We note that A/"m(£) decreases monotonically with increasing £ > 0, and < A/" m (£) < 1 for 
m > 1; for m — 0, A/i)(f) is positive and can be greater than 1 when £ becomes sufficiently small. 

The choice of unaligned perturbations is not unique, and the perturbation potential-density pairs may involve j3 parameter 
such that the background plus coplanar MHD perturbations are altogether scale-free (e.g., Syer & Tremaine 1996). As an 
example of illustration, we here presume coplanar logarithmic spiral perturbations defined by the density-potential pair l|55fl 
in the following analysis. For the specific /3 = 1/4 case, the background surface mass density and the perturbed surface mass 
density bear the same radial scaling oc r~ 3//2 . 



3.2.1 Marginal Axisymmetric Stability Curves 

We start from the axisymmetric coplanar MHD perturbation equation 140H with lu 7^ 0. For radial oscillations with a logarith- 
mic potential-density pair, we presume a surface mass density S oc r _3 / 2 + 1 ? ; a gravitational potential V = —2-KGMorS and 
thus a radial speed U oc r 1 / 2 + 1 f+ 2 according to the first of equation 14011 . After proper rearrangements, the set of equations 
14011 can be combined to form a single equation in terms of U(r), namely 



(UJ 2 - K 2 )U 





a 2 


[( 




a 2 






{[ 



4 J r 2 \ 4 



u 

(57) 



2+4/3- 

where C = (1 + 2/3)/ (2(3To), a 2 /r 2 = k 2 /[(1 + 2y9)(2 - 2/3)Z> 2 ] and lu 2 is extremely small. 

By inspection, equation l|57ft bears similar features of the dispersion relation in the classic WKBJ regime involving 
magnetic field. As the right-hand side of equation Ij57^ is real, axisymmetric instabilities first set via neutral lu = modes. It 
should be emphasized that the scale-free condition is met only for neutral modes. With to 2 = 0, the marginal stability curves 
are then given by 

D 2 = {1 - FCNq[1 - (1 - 4/% 2 /(2 + 4/3)]}(£ 2 + 1/4) + q 2 (e - P + 1/4) 
FCA/- (£ 2 + l/4) + (l + 2/3)(-2 + 2/3) ' 

As a check of necessary consistency, the axisymmetric marginal stability curve for an MSID when j3 = (C = 1) then gives 

n2 _ [l + g 2 -FA/- (l-g 2 /2)](g 2 + l/4) 
FA/- (£ 2 + l/4)-2 ■ 

One can readily show* that stationary dispersion relations 1581 and 1591 automatically satisfy inequality 1171 for physical 
solutions D 2 > 0. In the MSID case with different values of q parameter for magnetic field strengths, one readily obtains 
the marginal stability curves that separate two distinct unstable regimes, namely, the collapse regime for large radial spatial 
scales and the ring fragmentation regime for relatively short radial wavelengths (Lemos et al. 1991; Shen & Lou 2003, 
2004a). There exists a critical £ = £ c where D 2 diverges and this £ c is determined by FW*o(£ 2 + 1/4) = 2. For ( < we 
have -F7Vo(£ 2 + 1/4)— 2<0 and thus an enhancement of q tends to suppress the collapse boundary. For ( > we have 
FAfo (£ 2 + 1/4) — 2 > and thus an enhancement oi q tends to raise the ring fragmentation boundary. In short, an enhancement 
of ring magnetic field tends to reduce dangers of both collapse and ring fragmentation instabilities in the MSID case. 

While the enhancement of ring magnetic field tends to suppress ring fragmentation instabilities in general, it only sup- 
presses collapse instabilities for — l/4</3<l/4 and tends to aggravate collapse instabilities for 1/4 < /3 < 1/2 (see Fig. 2 for 
details). It has been shown that the magnetic force is radially inward and outward for —1/4 < j3 < 1/4 and 1/4 < P < 1/2, 
respectively. Why does the trend of variations for the collapse-stability appear seemingly paradoxical? This situation can 
be understood because the unperturbed background is in an MHD radial force balance and the surface mass density, the 
rotation speed and the magnetic field are coupled explicitly through condition llbH . Furthermore, for fast MHD density waves 
with m = (Lou & Fan 1998a), the gas pressure and magnetic pressure together is associated with the radial wavenumber 
squared, while the surface mass density and self-gravity are associated linearly with the radial wavenumber (see the first line 
of dispersion relat ion 1571 with small w 2 ). For —1/4 < f3 < 1/4, the surface mass density becomes smaller for stronger magnetic 
field strength, while for 1/4 < f3 < 1/2, the situation reverses. For the ring fragmentation instability at relatively large radial 
wavenumbers, the dominant magnetic pressure effect tends to stabilize the disc. For the Jeans collapse instability at relatively 
small radial wavenumbers, the dominant self-gravity effect is proportional to the background surface mass density; the roles 
of magnetic field for two different situations can then be readily understood. It is the coupling effect between the surface mass 
density and the magnetic field of the background that gives rise to this collapse feature. 

t We recall that inequality 1171 is automatically satisfied for l/4</3<l/2 and we need only to examine the range of— l/4</3<l/4. 
When the denominator of the right-hand side of equation 1581 is negative, the numerator must also be negative for D 2 > 0. It then 
follows that 1 — (1 — 4/3)g 2 /(2 + 4/3) > and inequality 1171 is met. When the denominator of the right-hand side of equation 1581 is 
positive, it is easy to show that inequality 1171 holds. As solution 1591 is a special case of solution 1581 . inequality 1171 is thus met. 
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The marginal stability curves for scale- free discs with /3 S ( — 1/4, 1/2) and without magnetic fields are 
(l-FCA/'o)(£ 2 + l/4) 



(60) 



FCAfoiZ 2 + 1/4) + 2(2/3 + l)(/3 - 1) ' 
that are consistent with those for the single-disc case of Shen & Lou (2004b; e.g., see their subsection 3.2.1). 

The marginal axisymmetric stability curves generally consist two branches, namely, the collapse branch and the ring- 
fragmentation branch (Lemos et al. 1991; Shu et al. 2000; Lou 2002; Lou & Shen 2003; Shen & Lou 2003). The lowest value 
of D 2 for stability is the maximum of the collapse branch and the highest value of D 2 for stability is the minimum of the 
ring- fragmentation branch. We note that the maximum of the collapse branch is always located at £ = for > —0.104, while 
for f3 < —0.104 the maximum of the collapse branch may locate at £ > with £ = corresponding to a local minimum of D 2 . 
It was proven in Appendix C of Shen & Lou (2004b) that £ = is always a local extremum for D 2 . Variations of the stable 
range for D 2 with parameters F, /3 and q are shown in panels (a) and (b) of Fig. 2. Finally we note that the axisymmetric 
destabilization involves only the fast MHD wave branch because the slow MHD wave branch is negative and thus unphysical. 
For non-axisymmetric (m 7^ 0) perturbations, both fast and slow MHD density waves are possible and stationary global 
configurations may represent transitions from stable to unstable situations (e.g., Shu et al. 2000). Here, the stability of our 
slow MHD density waves differs from those models involving finite disc thickness where either vertical or horizontal weak 
magnetic fields may induce Velikhov-Chandrasekhar-Balbus-Hawley instability and magneto-rotational instabilities (MRI) 
through slow MHD modes (Chandrasekhar 1961; Balbus & Hawley 1998; Kim & Ostriker 2000; Lou, Yuan & Fan 2001). 

It would be interesting to recall the axisymmetric stability analysis in the WKBJ approximation. The local dispersion 
relation for fast MHD waves in an azimuthal magnetic field is (Fan & Lou 1996) 

(u> - mfl) 2 = k + k 2 {C 2 A + a 2 - 27rGE /|fc|) , (61) 

where k is the radial wavenumber and other notations have the conventional meanings. Equation 10 U gives the marginal 
stability curve for axisymmetric m = perturbations as 

k 2 + k 2 (C 2 A + a 2 - 27rGE /|fc|) = . (62) 

By using scale-free disc equilibrium condition l|16|l and inserting £ = \k\r for large k and £ into dispersion relation i10'2fi . we 
readily obtain the corresponding marginal stability curve in the WKBJ regime as 

n 2 = (1 + g 2 )e - FC[1 - (1 - 4/3)g 2 /(2 + 4/3)]g 

FC£ + (l + 2/3)(-2 + 2/3) ' { °°> 

which bears a strong resemblance to expression l|58t of our global analysis, especially in view of the asymptotic expression for 
Afm(£,), namely 

A/- m (£) « (m 2 +£ 2 + l/4)- 1/2 (64) 

when m 2 + £ 2 2> 1. Another equivalent way of this comparison is to note the parallel terms in the WKBJ condition 1)62^ and 
in the first line of the corresponding global condition 1571 . also in reference to the above asymptotic expression for A/" m (C)- 
A specific illustrating example of such a comparison is shown in Fig. 3 with parameters /3 = (C = 1; flat rotation curve), 
F = 1 (full disc) and q = 0.1. In the short-wavelength limit (i.e., £ 3> 1), the WKBJ approximation is well justified and the 
results are qualitatively consistent, while in the long-wavelength of Jeans collapse regime, the WKBJ analysis differs from the 
global exact treatment significantly. 



3.2.2 Unaligned Logarithmic Spiral Cases with m > 1 

In parallel with the aligned m > 1 cases, we now construct unaligned logarithmic spiral solutions from coplanar MHD 
perturbation equations (1361) — 1391) . For the density-potential pair S oc r~ 3//2+1 ^ and V = —2%GJ\f m rS, we then have radial flow 
speed U oc r _1 / 2 +' 3 + 1 5 and z— component specific angular momentum J oc r 1 / 2 +' 3 + 1 ?. Coplanar MHD perturbation equations 
lp |l - P ty can therefore be written in the forms of 

milS H ■ 2joh7 H r — J — , 



20. 

mQ.iU H J = 

r 



r£ 



- 2ttGA/; 



(1 - 4(3)C 2 A 



mQJ + r (1 - p)mU = m { 2nGMm - — )rS ■ 
V rZ J 



2E r 
Cl(l-4/3) 
2Q,r 



S- 



■mflr 2 



iU 



B iU 



R = 



B U 



mQ,r fir 
The first three equations of set 165H immediately lead to 

aiS + bull + ciJ = , a 2 S + b 2 iU + c 2 J = , 



2-2/3)i£-r 



a 3 S + b s iU + c 3 J = 



iU 



(65) 



(66) 



© 0000 RAS, MNRAS 000, 000-000 



12 Y. Shen, X. Liu, Y.-Q. Lou 




Figure 2. Variations of the stable D 2 range for (a) full discs with F = 1 and (b) f3 = 1/4 discs. The solid, dashed and dotted curves 
stand for q = 0, 0.3 and 1, respectively. For specific ranges of either (3 £ (—1/4,1/2) [see panel (a)] or F g (0,1] [see panel (b)], a 
horizontal line intercepts the stability boundary curves at two values of D 2 and the stable range of D 2 lies between the two intersection 
points. The D 2 range to the left of the smaller intersection point corresponds to the collapse regime and the D 2 range to the right 
of the larger intersection point corresponds to the ring-fragmentation regime. Panel (a): For full discs with F = 1, an enhancement 
of magnetic field tends to suppress ring-fragmentation instabilities in general. Meanwhile, an enhancement of magnetic field tends to 
suppress collapse instabilities for —1/4 < f3 < 1/4 but to aggravate collapse instabilities for 1/4 < fi < 1/2 and bears no effects on 
collapse instabilities when (3 = 1/4. For stability curves of moderate magnetic field strengths (e.g., the solid and dashed curves), the 
disc system becomes less stable for both instabilities when (3 decreases. As (3 becomes sufficiently small, there will be no stable range 
for D 2 [e.g., for q = without magnetic field, this lower limit of fi is ~ —0.130 as derived in Appendix D of Shen & Lou (2004b)]. For 
a stronger magnetic field (e.g., the dotted curve), the trend for ring-fragmentation instabilities remains the same while the trend for 
collapse instabilities is different depending on whether (3 < 1/4 or not; the collapse regime may disappear completely when (3 becomes 
small enough. Panel (b): For f3 = 1/4 discs, an enhancement of magnetic field tends to suppress the ring-fragmentation regime but does 
not influence the collapse regime because the maximum of the collapse branch remains always at £ = such that the heavy solid curve 
on the left actually represents a superposition of the three q curves. A decrease of F (i.e., a more massive dark-matter halo) tends to 
suppress both instabilities. 
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Figure 3. A comparison of the global axisymmetric stability treatment (solid curves) with the local WKBJ approach (dashed curves) 
for the case of f3 = 0, F = 1 and q = 0.1. The marginal stability curves for the long-wavelength Jeans collapse and the short-wavelength 
ring fragmentation regimes arc denoted. For £ 3> 1 (i.e., the short-wavelength limit), the WKBJ approximation is well applicable, while 
for small £ the results in the WKBJ approximation deviates from those of the exact global treatment considerably (see Shen & Lou 2003, 
2004b). 
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where coefficients ai, bi and Cj (i — 1, 2, 3) are explicitly defined by 



ai = mO, , 

_ (1 - 40)C 2 A 
a 2 - ^ 

a 3 = — mr[ 2-kGM„. 



oi = 



mEn 
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27rGMr, 



Cl 
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6 3 = r(l-/3)n- 
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(l-4/3)Cj 
2«r 



[(-2/3 + l/2)ig - £ 2 - m 2 ]Cj 
mQ.r 2 

, c 3 = mfi . 



C2 



2Q 
r 



(67) 



For nontrivial solutions of the set of homogeneous algebraic equations 1661 , the determinant must vanish 

bi(a 2 c 3 - a 3 c 2 ) c 1 (a 2 b 3 - a 3 b 2 ) ,„ Q . 

ai 1 1 1 — I 1 = u • 

O2C3 — b 3 c 2 b 2 c 3 - b 3 c 2 

By using the set of expressions 1671 and applying equilibrium conditions 1161 . we obtain the final dispersion relation for 
stationary logarithmic spiral m > 1 MHD configurations in the form of a quadratic equation in terms of y = D 2 , namely 

C 2 y 2 + Ciy + Co = , (69) 

where the coefficients C 2 , C\ and Co are explicitly defined by 

C 2 = Hm , 

"(3-4/3)m 2 -8/3 2 - (4£ 2 - 3)/3 + 3^ 2 - 1/4 ^„. r . 2 . .2 7 



<h = - 



C 



2(1 + 2/3) 



FCN m -- + 73 



q - Am (I — FCNrr, 



(70) 



(1-4/3) 



f-2m 2 



4(1 + 2/3) 

with auxiliary parameters defined by 
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-2g 2 - 2/3 + 1/2) ' 

(TT^fl) ™ m + 1 " 4/3 . 



4 _ (-2m 2 - 2g 2 - 2/3 + 1/2) 



An = m 



- f + \ + 2/3 



B m = (1 + 2/3) (m -2 + 2/3) 
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2(1 + 2/3) 

(1 + 2/3) 
2/3P 



(1 - FCM m )q 



7~tm — FCJ\f m ^A m -\- J3n 



(71) 



that are the same as those in Shen & Lou (2004b) for unaligned logarithmic spiral cases. For a positive determinant A 
C\ — 4C2C0 > 0, the two real solutions to quadratic equation 169H are given by 



3/i = 



-d+A 1 / 2 



and 



-d - A 1 / 2 



■ic, " —Ic—- ^ 

Similar to m > 1 cases for aligned perturbations in the limit of vanishing magnetic field (i.e., q — » 0), equation 1691 
reduces to 

H m y 2 - A m (l - FCAf m )y = 
that leads to a nontrivial solution 

n2 (1 - FCAf m )Am 

D = 2/1 = - 



Tin 



(73) 



(74) 



Solution 17411 is exactly the same result for the single-disc case discussed in Section 3.2 of Shen & Lou (2004b). Another 
limiting regime of — > has been analyzed by Lou (2002) for global stationary perturbation structures in a single MSID. 



3.2.3 Force-Free Magnetized Discs with (3 = 1/4 

For stationary spiral MHD configurations, we here analyze again the = 1/4 case as a specific example of illustration. The 
background equilibrium ring magnetic field is force-free and inequality 1171 is satisfied. The determinant A of quadratic 
equation 1691 remains always positive as shown in Appendix A and hence two real solutions of y — D 2 to equation 1691 exist. 
In parallel with the aligned = 1/4 case, we solve quadratic equation 1691 for given values of parameter F and display in 
diagrams the two branches of solutions j/i and y 2 (y = D 2 ) versus £ variation. 

As Si < for m = 1 and when F becomes sufficiently small, there is a diverging point for y\ = D 2 solution where 
C 2 — Hi = 0. Since Tii increases monotonically with increasing £ for (, > and attains its minimum value at £ = (Hi is 
an even function of £; see Appendix C of Shen & Lou 2004b), the existence of such a critical point £ c for a diverging y = D 2 
requires 

FCJVi(0){5/4 + 2/3) + 40 2 - 1 < . (75) 

For the specific — 1/4 case, this inequality 1751 implies a critical F value of F c = 3/[7CA/i(0)] below which there exists 
one divergent point of £ c where yi diverges. Numerically, we have determined F c — 0.6842 [note that this F c is of the same 
value as that of the aligned m = 1 case because A/"i(0) = 7 3 i(l/4)]. For selected F values in the range of < F < F c , the 
corresponding critical £ c at which TLi = can be numerically computed; for example, we have £ c = 1.428 for F = 0.5. 

It becomes much simpler for cases of m > 2 because C 2 remains always positive and there is no diverging point for D 2 
solutions. Therefore, yi and y 2 remain always as the upper and lower branches, respectively. For m = 1 and m — 2 cases with 
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Figure 4. Two D 2 solution branches yi and j/2 of quadratic equation 1691 for m = 1, 2, P = 1/4, F = 1, 0.5 and 5 = 0, 1. The y\ and 
j/2 branches are plotted in solid and dashed lines, respectively. Parameters adopted are shown in panels (a) — (c). 



different parameters of F and q, we show in four panels (a) — (d) of Fig. 4 relevant curves of y = D 2 versus £ variation to 
illustrate the basic features. To be more explicit for the case of /3 = 1/4 in quadratic equation 1691 . we write 



FCNm 



.2 3 

■^+4 



2 lm 



2_ 3 
2 



+ (m 2 + f + ^ (1 - FCNm) 



D 2 + -(m 2 + - FCAf m )q 2 = , 
where auxiliary parameters Af m , Am, B m , C and 7i m are all evaluated for the special value of (3 — 1/4. 



(76) 



4 PERTURBATION PHASE RELATIONSHIPS AND ANGULAR MOMENTUM TRANSFER 

For multi- wavelength observations of nearby spiral galaxies (Sofue et al. 1986; Kronberg 1994; Beck et al. 1996), it is possible 
to identify large-scale spatial phase relationships of spiral patterns using various observational diagnostics (e.g., Mathewson 
et al. 1972; Visser 1980a, b; Neininger 1992; Beck & Hoernes 1996; Frick et al. 2000, 2001; Ferriere 2001; Lou et al. 2002). For 
this reason, we derive below spatial phase relationships for various coplanar MHD perturbation variables, keeping in mind of 
potential applications to magnetized spiral galaxies (e.g., Fan & Lou 1996; Lou & Fan 1998a). 
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4.1 General Perturbation Phase Relationships for both Aligned and Spiral Cases 



Here we systematically examine the spatial phase relationships among enhancements of gas surface mass density, magnetic field 
and velocity perturbations by analysing the relationships among S, Z, R, J and iU for aligned and logarithmic spiral cases, 
separately. Starting from the set of homogeneous algebraic equations 1 1H for the aligned cases, we obtain the relationships 
among S, U and J in the following, namely 
iU C3CL1 — C1C13 J 0261 — ai&2 

~S = 



C1&3 - C3&1 
where at, bi and Ci (i 

iU mQr . 
~~S 



S C162 — C261 

1, 2, 3) are defined by expressions 14511 . By using definitions 1451 and equations 143H . we obtain 



(77) 



-0A 



z 
s 



(1/2 - 2/3)B 



©A , 



iR 
S 



"a, -x = — =— (1 - W<oa) 

S 2^0 



So S So S So 

for the aligned cases, where the common dimensionless real factor Qa is defined below by 

a 2 - 27rGP m E r - Q 2 r 2 



9 



(l/2-2/3)C2 -(l + 2/?)tt 2 r 2 



(78) 



(79) 



Similarly, by using the set of homogeneous algebraic equations 1661 . definitions 16 11 and equations 165H . we readily derive 
for the cases of global stationary logarithmic spiral configurations 



S 



mQr 



Z 

s 



iCBo 



iR 
S 



mBc 



e, 



So S So S So o 

where the common dimensionless complex factor Q$ is defined below by 
a 2 ~ 27rGW m E r - ftV 



Q.T 2 
' So 



e 



(80) 



(81) 



(1/2 - 2f3)C\ - (3 - 2i£)n 2 r 2 ' 
The two sets of algebraic equations 17811 and 18011 are the most general phase relationships for aligned and logarithmic 
spiral cases, respectively. Once a j3 value is specified and a particular D 2 for a stationary solution inserted, we immediately 
obtain the phase relationship between any two perturbation variables. As an example of illustration, we examine again the 
special /3 = 1/4 case and determine the relevant phase relationships in the next subsection. 



4.2 The f3 = 1/4 Force- Free Case 

When j3 = 1/4, the results of phase relations derived in the last subsection 4.1 can be simplified considerably as all terms 
involving the factor 1/2 — 2/3 vanish. Meanwhile, inequality l|17|l is automatically satisfied. The only requirement for a physical 
D 2 solution is simply D 2 > 0. As already noted in subsections 3.1.3 and 3.2.3, there are two branches of D 2 solutions. We insert 
relevant physical D 2 solutions in equations 1781 and 1801 to determine phase relationships among perturbation variables. In 
particular, we are interested in the phase relationship between the perturbed surface mass density and the perturbed azimuthal 
magnetic field. We now analyze below this mass-magnetic field phase relationships for aligned and spiral cases, separately. 
Some illustrative examples are shown in Fig. 5 and Fig. 6 to provide more direct visual impressions. 



4-2.1 Aligned Cases 

For aligned (3 = 1/4 cases, we readily find Z = which means the azimuthal magnetic field remains unchanged. From equation 
17811 . we derive 

iR/S oc 3D 2 /2 - 1 + FCV m {l + D 2 ) = (FCT m + 3/2)D 2 - (1 - FCT m ) . (82) 
Meanwhile, equation 1541 can be rewritten in the form of 

m 2 g 2 {l + 2FCP m /3) = D 2 {D 2 - (1 - FCV m )/[FCV m + 3(m 2 - 3/2)/(2m 2 + 3/2)]} 

(m 2 + 3/4:)[FCVm + 3(m 2 - 3/2)/(2m 2 +3/2)] D 2 - (1 - FCV m )/(FCP m + 3/2) ' ( ' 

that contains two distinct D 2 solutions according to subsection 3.1.3, namely, the plus-sign solution yi and the minus-sign 
solution 7/2 in the solution form I51H . respectively. As shown in Appendix B, we emphasize that the plus-sign solution y\ (if 
positive and thus physical) makes D 2 — (1 — FCV m ) / (FCV m + 3/2) > and the minus-sign solution 1/2 makes D 2 — (1 — 
FCPm) I {FCVm + 3/2) < on the right-hand side of equation 1831 1. respectively. Therefore, the two branches of D 2 solutions 
(if positive and thus physical) will give either iR/S > for the plus-sign j/i solution or iR/S < for the minus-sign i/2 solution. 

4-2.2 Cases of Stationary Logarithmic Spiral Configurations 

For spiral cases with [3 — 1/4, we readily obtain from equations I80H and 1811 that 

9s oc [(FCN m + 3/2)D 2 - (1 - TCM m )](3 + 2if) . (84) 
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(a) Fractional perturbed surface mass 
density gray-scale plot as a spatial 
phase reference for either fast or slow 
MHD density wave modes 
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(b) Fractional perturbed radial mag- 
netic field in a gray-scale plot for the 
fast MHD density wave mode 
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(c) Fractional perturbed radial mag- 
netic field in a gray-scale plot for the 
slow MHD slow density wave mode 



Figure 5. Spatial phase relationship between the perturbed surface mass density and the radial magnetic field for the aligned m = 2 
bar-like case for the specific case of /3 = 1/4. Shown in the figures above are the fractional perturbed surface mass density and magnetic 
field contours (i.e., perturbation variables divided by the radial scalings of their equilibrium quantities with r~ 3 / 2 for the surface mass 
density and r~ 1 for the azimuthal magnetic field). The brightest regions correspond to the strongest density enhancement. Panel (a) is 
for the perturbed relative surface mass density; panels (b) and (c) are for the perturbed relative radial magnetic fields for the fast and 
slow MHD perturbations respectively. According to equation 1821 . there is a phase shift of ^tt/2 for the global stationary fast and slow 
MHD density wave modes. Since m = 2, the pattern appears a rotation of ^fir/A. 

Following the same procedure, we can show that the two branches of D 2 solutions, namely, the plus-sign solution j/i and the 
minus-sign solution y-2 given by equation (1721 . if positive and thus physical, correspond to the factor (FCAfm + 3/2) D 2 — 
(1 — FCAfm) being positive and negative, respectively (see Appendix B for a detailed proof). We now examine the phase 
relationship between coplanar perturbation enhancements of surface mass density and magnetic field. 

In the short- wavelength limit (i.e., ( > 1 with a logarithmic spiral pattern tightly wound), we have Z/S oc ±1 and 
R/S oc ±1 where the plus- and minus-signs correspond to the y\ and j/2 solutions, respectively. Therefore, the phase relations 
between perturbation enhancements of surface mass density and magnetic field are in-phase for the yi solution (if positive 
and physical) and out-of-phase for the t/2 solution. For m > 2, yi remains always greater than j/2- In contexts of magnetized 
spiral galaxies, these solution properties in the WKBJ regime are consistent with the physical identifications of the y\ mode 
with a stationary fast MHD density wave and the 1/2 mode with a stationary slow MHD density wave, respectively (Fan & 
Lou 1996; Lou & Fan 1998a; Lou & Fan 2002). It is now possible to model magnetized spiral galaxies in terms of coplanar 
stationary fast and slow MHD density waves with more general rotation curves including the case of a flat rotation curve for 
a magnetized singular isothermal disc (Shu et al. 2000; Lou 2002; Lou & Fan 2002; Lou & Zou 2004a, b; Lou & Wu 2004). 

In the long-wavelength limit (i.e., f < 1 with a logarithmic spiral pattern widely open), we have Z/S oc =Fi and R/S oc =pi 
where the minus- and plus-signs correspond to the j/i and 2/2 solutions, respectively. Therefore, the coplanar perturbation 
enhancement of magnetic field is either ahead of or lagging behind the perturbation enhancement of surface mass density by 
a phase difference of ~ 7r/2. 

For intermediate radial wavelengths, the phase difference between perturbation enhancements of surface mass density 
and magnetic field in a logarithmic spiral pattern can be readily determined for given specific conditions. 

On the basis of our model analysis here, it appears that except for global stationary fast MHD density waves in the 
extreme WKBJ regime, magnetic and gas spiral arms tend to be phase shifted relative to each other for a global stationary 
logarithmic spiral pattern for a much broader set of rotation curves including the special case of a flat rotation curve. 

4.3 Angular Momentum Conservation 

One can readily show that for aligned global stationary MHD perturbation configurations, there is no net radial angular 
momentum flux as aligned stationary MHD density waves do not propagate radially (Shu et al. 2000; Lou 2002). We here 
briefly discuss the process of angular momentum transfer in association with the stationary logarithmic spiral MHD modes 
obtained in this paper. According to the analyses of Lynden-Bell & Kalnajs (1972) and Fan & Lou (1999), the total angular 
momentum flux carried by coplanar MHD density waves should contain three separate contributions, namely, the advective 
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(a) Fractional perturbed mass density 
in a gray-scale plot as a spatial phase 
reference for either fast or slow MHD 
density wave modes 



(b) Fractional perturbed azimuthal 
magnetic field in a gray-scale plot for 
the fast MHD density wave mode 



(c) Fractional perturbed azimuthal 
magnetic field in a gray-scale plot for 
the slow MHD density wave mode 



Figure 6. Phase relationship between the perturbed surface mass density and the azimuthal magnetic field for the logarithmic spiral 
m = 2 case (an effective radial wavenumber § = 5) for the specific case of = 1/4. Shown in above three panels in order are the 
fractional perturbed surface mass density and azimuthal magnetic field gray-scale plots (i.e., divided by the radial scaling of their 
respective equilibrium quantities with r~ 3 ' 2 for the surface mass density and r~ 1 for the azimuthal magnetic field). The brightest 
regions correspond to the most density enhancement. Panel (a) is for the perturbed relative surface mass density; panels (b) and (c) 
are for the perturbed azimuthal magnetic fields for the fast and slow MHD perturbations, respectively. According to equations 1841 and 
ISTfl . there is a phase shift of ~ — 0.0937T (~ 0.9077r) for the fast (slow) MHD perturbation mode. Since m = 2, the pattern appears a 
rotation of half of the phase shift. For the fast MHD density wave mode, the magnetic field enhancement roughly follows the spirals of 
surface mass enhancement. For the slow MHD density wave mode, the density spiral arms and the magnetic spiral arms largely interlace 
with each other. For f > 1, the phase pattern relationships become those in the usual WKBJ regime (e.g., Fan & Lou 1996; Lou & Fan 
1998a). 
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The perturbed three-dimensional gravitational potential <j>\(r,6,z) associated with a Fourier component of a coplanar 



logarithmic spiral perturbation in surface mass density Ei = err 
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is then 



fa{r,B,z) = ~2ixGoe- Lma I dfce" fc|z| J m (kr) r" 1/2+1 « J m (kr a )dr a , 
Jo Jo 

where J m (u) is the cylindrical Bessel function of order m with an argument u (e.g., Binney & Tremaine 1987). In fact, integral 
l|88|l gives a simple form of Kalnajs function in the disc plane at z — with the real part of <f>\{r, 8, z) explicitly given by 



i) = -2nGo / dfce" fc|z| J m (fer) 



-1/2 



cos(m# - (,lnr a )J m (kr a )dr a 



(89) 



io Jo 

It is fairly straightforward to show that A A and A s are independent of r where the phase relationships for stationary 
logarithmic spiral cases derived in Section 4.1 have been use. While mathematically tedious, one can further show that A G 
is also independent of r (see Appendix C). Therefore, the total angular momentum flux A' 4 + A G + A B remains constant at 
all radii, a natural result of angular momentum conservation given that perturbation surface mass density falls off as r~ 3 / 2 . 
Therefore, in order to maintain a stationary logarithmic spiral pattern, there must be a source or a sink at the disc center. The 
perturbation energy flux, however, is zero at all radii for the stationary coplanar MHD configurations [see equation (4.17) of 
Fan & Lou (1999) in the WKBJ regime]. The mass flux is also zero at all radii since the first-order mass flux J Q * EoMirdf? = 0, 
and the second-order mass flux f£* (Yj-lui +EoU2)rd0 = [where U2 is the second-order perturbed radial velocity, this equation 
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holds due to the second-order mass conservation (Fan & Lou 1999)]. Therefore our model result differs from that of Spruit 
(1987) for stationary logarithmic spiral shocks, where only the gravity of the central object is considered. In Spruit's model, 
there is a net mass accretion through the disc and an effective viscosity parameter "a" (Shakura & Sunyaev 1973) may be 
estimated. 

Physically, the central source of angular momentum required to sustain global stationary logarithmic spiral patterns may 
be plausibly identified with a fast spinning supermassive black hole (SMBH) in the nucleus of a spiral galaxy. The Kerr 
spacetime associated with such a spinning SMBH excites and sustains MHD density waves in the surrounding disc through a 
long-range gravitational interaction. For trailing spiral patterns, a net flux of angular momentum is transport radially outward 
in the disc. Conceptually, this process should cause the spinning SMBH to gradually slow down, entirely analogous to the 
spin-down of a rotating neutron star by emitting electromagnetic waves that carry a net angular momentum away. 



5 SUMMARY 

In this paper, we have explored and analyzed stationary perturbation structures in coplanarly magnetized razor-thin scale- 
free discs in a more general manner. To be specific, we start from a rotational and magnetized background equilibrium of 
axisymmetry that is dynamically self-consistent with a surface mass density Eo oc r~ a , a rotation curve vq oc r _/3 , a purely 
azimuthal (ring) magnetic field Bo oc r -7 with a vertically integrated barotropic equation of state in the form of n = KT. n . 
The radial force balance at all radii in an MHD disc implies several simple relationships among these power indices a, f3, 7 and 
n, as stated explicitly in equation l)15|l . Without loss of generality, we can simply use f3 as an independent parameter to specify 
properties of an important class of radial variations for rotational MHD background equilibria. The allowed range of f3 falls 
within the interval (—1/4, 1/2) (Syer & Tremaine 1996; Shen & Lou 2004a, b) with the special case of /3 = corresponding to 
a magnetized singular isothermal disc (MSID) system (Shu et al. 2000; Lou 2002; Lou & Fan 2002; Lou & Zou 2004a, b; Lou & 
Wu 2004). For clarity and convenience of our analysis, we introduced several dimensionless parameters. The first parameter is 
the partial disc parameter F = <j>/<f>T for the ratio of the gravitational potential ((f)) arising from the disc to that (4>t) arising 
from the entire system (Syer & Tremaine 1996; Shu et al. 2000; Lou 2002) including an axisymmetric halo mass distribution 
that is involved in the background equilibrium but is presumed to be unresponsive to coplanar MHD disc disturbances (a 
massive dark matter halo is an example in mind); for a partial disc system, we have < F < 1, while for a full disc system, we 
have F = 1 by conventional definitions. The second parameter is an effective rotational Mach number D defined as the ratio 
of the disc rotation speed to the sound speed yet with an additional scaling factor (1 + 2/3) _1/ ' 2 ; D is constant at all radii. The 
third parameter q is a measure for the background ring magnetic field strength defined as the ratio of the azimuthal Alfven 
speed to the sound speed; and q remains constant at all radii. For 1/4 < (3 < 1/2, there is no constraint on q and the net 
Lorentz force of the background azimuthal magnetic field is radially outward (the magnetic pressure force is stronger than the 
magnetic tension force) to resist the self-gravitation; for (3 = 1/4, the background azimuthal magnetic field is force-free (the 
magnetic pressure and tension forces cancel each other exactly); for — 1/4 < /3 < 1/4, the net Lorentz force of the background 
azimuthal magnetic field points radially inward (the magnetic pressure force is weaker than the magnetic tension force) to 
aggravate the self-gravity and thus a restriction [i.e., inequality l|17fl ] is imposed on the magnetic field strength in order to 
guarantee that the surface mass density be positive in equation 11611 . Inequality 1171 provides a necessary physical criterion 
for plausible D 2 solutions. 

With such a rotational MHD background equilibrium chosen, we introduce coplanar MHD disturbances to construct 
global stationary perturbation structures as viewed in an inertial frame of reference. Perturbation variables are expressed in 
terms of Fourier components for either aligned or logarithmic spiral pattern forms (Syer & Tremaine 1996; Shu et al. 2000; 
Lou 2002; Lou & Shen 2003; Shen & Lou 2003; Lou & Zou 2004a). We derive analytical solutions for stationary coplanar 
MHD perturbation configurations in both aligned and logarithmic spiral cases. We now summarize key results below. 

(i) Cases of Aligned Coplanar MHD Perturbations 

For aligned cases, we choose perturbation variables that bear the same power-law dependence in radius r as the background 
MHD equilibrium does. For the special aligned case of m = 0, the stationary coplanar MHD perturbations actually describe 
alternative equilibria to the axisymmetric background equilibrium with proper rescaling factors (Shu et al. 2000; Lou 2002; 
Lou & Zou 2004a) and are therefore somewhat trivial. 

For aligned cases of m > 1, we derive a quadratic equation 1481 in terms of y = D 2 . The resulting positive D 2 is the rotational 
parameter required for sustaining a stationary nonaxisymmetric MHD perturbation configuration. Other parameters involved 
in quadratic equation l|48[l include m, F and q. For a general combination of /3, m, F and q, the determinant A of quadratic 
equation l|48[l may not be always positive definite. On the other hand, for specific values of /3, equation Ij48^ always gives two 
different real solutions (see Appendix A). Once we prescribe parameters /3, m, F and q, the two solutions of D 2 can be readily 
determined. For physical solutions, D 2 must be positive and satisfy inequality (1171 simultaneously. 

As a specific case study, we examined the case of (5 = 1/4. Here the background azimuthal magnetic field remains force-free 
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and the D 2 solution is free from constraint (I17H . It has been further shown in Appendix A that for this f3 = 1/4 case, equation 
14811 always has two distinct real solutions, the plus-signed solution y\ and the minus-signed solution yi, as contained in 
equation 15H . All these aspects contribute to a significant simplification of the relevant computational procedures. We derive 
two branches of D 2 solutions and need to require only D 2 > 0. There is yet a qualitative difference between the m = 1 case and 
the m > 2 cases. For m — 1, there exists one critical value F c — 0.6842; for F below F c , yi and yi are negative and positive, 
respectively, while F above F c , both branches of D 2 solutions are positive and thus physical. For m > 2 in comparison, both 
branches of D 2 solutions remain always positive and thus physical. Moreover, yi remains always to be the upper branch for 
azimuthal stationary fast MHD density waves and represents the hydrodynamic counterpart of the unmagnetized single disc 
case; complementarily, yi remains always to be the lower branch for azimuthal stationary slow MHD density waves (Fan & 
Lou 1996; Lou & Fan 1998a) and is caused by the very presence of the azimuthal magnetic field (detailed computations reveal 
that 2/2 must be smaller than 2/3 which seems to be too small for galactic applications). When a disc rotates sufficiently 
fast, it could only support the upper y\ stationary MHD configurations. An increase of magnetic field parameter q will raise 
both physical branches of D 2 solutions as displayed in Fig. 1. The variation trend of the upper y\ branch appears to be more 
sensitive. Therefore for a much stronger magnetic field, a sufficiently rapidly rotating disc may sustain stationary coplanar 
MHD configurations in terms of the upper y\ = D 2 branch. 

(ii) Cases of Logarithmic Spiral Configurations for Coplanar MHD Perturbations 

For cases of stationary spiral perturbation structures, we choose coplanar MHD perturbations in terms of Kalnajs logarithmic 
spirals (Kalnajs 1971; Lemos et al. 1991; Shu et al. 2000; Lou 2002; Lou & Fan 2002; Lou & Shen 2003; Lou & Zou 2004a; Lou 
& Wu 2004). For the special m — case with radial propagations, it happens that inequality I17H is automatically satisfied 
and the stationary coplanar MHD configurations with y = D 2 > obtained here represent marginal stability curves (Syer & 
Tremaine 1996; Shu et al. 2000; Shen & Lou 2003, 2004b; Lou & Zou 2004b). Only when the rotational parameter D 2 falls 
within a specific finite range can the magnetized disc be stable against axisymmetric disturbances at all wavelengths. A disc 
with too slow a rotation speed will succumb to Jeans' instability in the collapse regime corresponding to long wavelength 
perturbations, while a disc with too fast a rotation speed will suffer the ring-fragmentation instability corresponding to 
relatively short wavelength perturbations (Safronov 1960; Toomre 1964; Lemos et al. 1991; Lou & Fan 1998a, 2000; Shu et 
al. 2000; Shen & Lou 2003, 2004a; Lou & Zou 2004b). In comparison to a hydrodynamics disc system, the enhancement of 
azimuthal magnetic field tends to suppress ring-fragmentation instabilities as expected from the perspective of the MHD Q 
parameter (Lou & Fan 1998a; Lou 2002; Lou & Zou 2004a, b), while an enhancement of azimuthal magnetic field tends to 
suppress Jeans' collapse instabilities for —1/4 < j3 < 1/4, to aggravate Jeans' collapse instabilities for 1/4 < (3 < 1/2 and 
to bear no effect on Jeans' collapse instabilities when f3 — 1/4. We have provided intuitive interpretations for the roles of 
magnetic field earlier. 

In parallel with the aligned cases, we derive quadratic equation l|69[l in terms of y = D 2 for logarithmic spiral cases of 
m > 1. In addition to the four parameters /3, m, F and q, we now have one more parameter £ for the radial wavenumber. 
While the determinant A of quadratic equation l)69|l may not be always positive in general, we prove in Appendix A that 
A > for the specific (3 = 1/4 case. Once the five parameters f3, m, F, q and £ are specified, the D 2 solutions can be readily 
obtained. For physical solutions, D 2 must be non- negative and also satisfy inequality 11711 . 

As an example of illustration, we study the specific /3 = 1/4 case inequality (1171 being satisfied automatically. Not surpris- 
ingly, the aligned cases and the logarithmic spiral cases parallel with each other very well, because of a mere change of wave 
propagation direction. For the special spiral case of m = 1, there again exists the same critical F c ~ 0.6842 for F value as 
in the aligned m — 1 case; for F < F c , there is one critical point of £ c where D 2 = y\ diverges. More specifically, we have 
unphysical yi < for < £ < £ c , while we have physical yi > for £ > £ c . There is no divergent point for D 2 solution 
branch and we always have yi > (yi = when q = 0). The solution structures become simpler for m > 2 cases when 
both D 2 branches of solutions are non-negative and there is no divergent point of D 2 = y±. Moreover, y\ remains always 
to be the upper branch for stationary spiral fast MHD density waves and is the hydrodynamic counterpart for those in the 
unmagnetized single disc case, while yi remains always to be the lower branch for stationary spiral slow MHD density waves 
(Fan & Lou 1996) and is caused by the very presence of the magnetic field (detailed computations reveal yi must be smaller 
than 2/3). 

(iii) Phase Relationships among Coplanar MHD Perturbation Variables 

So far, we have constructed global stationary coplanar MHD configurations for both aligned and logarithmic spiral cases 
in magnetized discs with a range of rotation curves. As expected, the inclusion of azimuthal magnetic field gives rise to one 
additional distinct physical solution for stationary slow MHD density waves to the problem (Fan & Lou 1996; Lou & Fan 
1998a; Lou 2002; Lou & Zou 2004a); the Alfvenic modes are excluded by our restricted consideration for coplanar MHD 
perturbations. The two D 2 branches of solutions are at least mathematically reasonable and may both have applications 
to astrophysical disc systems such as magnetized disc galaxies and so forth. These devised MHD disc problems are highly 
idealized with yet precious global analytical perturbation solutions, and conceptually still belong to a subclass of more general 



© 0000 RAS, MNRAS 000, 000-000 



20 Y. Shen, X. Liu, Y.-Q. Lou 



time- dependent MHD perturbation solutions. For applications to magnetized spiral galaxies, it would be of considerable 
interest to discuss pertinent physical aspects and identify relevant observational diagnostics. 

For physical solutions of D 2 > 0, the j/i branch would be applicable to MHD discs with relatively fast rotation while the 
2/2 branch would be applicable to MHD discs with relatively slow rotation in general. For spatial phase relationships, the 
two branches of D 2 > solutions give different results. We take the specific (3 = 1/4 case as an example of illustration. For 
corresponding aligned cases, the plus-solution y\ will lead to lR/S > and the minus-solution 3/2 will lead to iR/S < 0; 
there are no azimuthal magnetic field perturbations for aligned (3 = 1/4 cases. For tightly wound logarithmic spiral cases, 
perturbations of both radial and azimuthal magnetic field are approximately in phase with the perturbation enhancement of 
the surface mass density for the y\ branch, while they are approximately out of phase for the 1/2 branch. For open logarithmic 
spirals, the radial and azimuthal magnetic field perturbations are either ahead of or lagging behind the enhancement of the 
surface mass by a significant phase difference. These phase relationships together with those associated with flow perturbations 
(Visser 1980a, b) are valuable to interpret the spatial phase shifts between optical arms and magnetic arms in lopsided, barred 
or other spiral galaxies (Mathewson et al. 1972; Neininger 1992; Beck & Hoernes 1996; Beck et al. 1996; Fan & Lou 1996; Lou 
& Fan 1998, 2000, 2002; Frick et al. 2000, 2001; Lou et al. 2002). 

(iv) Angular Momentum Transfer in Steady Logarithmic Spirals for Coplanar MHD Perturbation Configurations 

Stationary logarithmic spiral configurations of MHD density waves carry constant angular momentum flux either outward 
or inward associated with the advective transport, the gravity torque and the magnetic torque (Lynden-Bell & Kalnajs 1972; 
Goldreich & Tremaine 1978; Fan & Lou 1999). Therefore, the net angular momentum flux in the entire disc system is conserved 
and there must be either a source or a sink at the disc center. For the y\ branch with D 2 > 0, we find that the total angular 
momentum flux is inward for leading (i.e., £ > 0) spiral MHD density waves and is outward for trailing (i.e., £ < 0) spiral MHD 
density waves. For the 1/2 branch in comparison, the total angular momentum flux is outward for leading (i.e., £ > 0) spiral 
MHD density waves and is inward for trailing (i.e., £ < 0) spiral MHD density waves. As a real MHD disc system typically 
rotates at a relatively large D 2 and inequality (1171 can also rule out small D 2 for (3 7^ 1/4, the lower 1/2 branch might be 
rarely useful in galactic applications. In either case, however, the net mass accretion rate is zero since if we focus on materials 
between two concentric circles, at any time the influent angular momentum from the inner ring is equal to the effluent angular 
momentum from the outer ring such that materials inside this belt neither gain nor lose any angular momentum and therefore 
no net mass accretion occurs. 
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APPENDIX A: BASIC PROPERTIES OF DETERMINANT A 



For m > 1, we first examine the solution criterion of quadratic equation 1481 for the aligned case. The determinant A of 
quadratic equation (1481 is 



A = CT - 4C 2 C 



(Al) 



where the three coefficients Cj (j — 0, 1, 2) are defined explicitly by 
"(3-4/3)m 2 - (1 - 4/3)(l - 6/3 + 4/3 2 ) 



C2 



2(1 + 2/3) 



-FCV m + mr-(l-4/3){2-j3) 



(1 ~ 4/3) 
(1 + 2/3) 



-2m 2 + (1 - 4/3) (1 - 2/3) 



(1 + 2/3) 



FCV m + 1-4/3 



d = 2(1- FCV m ){ A 



c = A 2 m {\-FCV, 



(3 - 4/3)m 2 - (1 - 4/3) (1 -6(3 + 4/3 2 



2(1 + 2/3) 
-2m 2 + (1 - 4/3) (1 - 2/3) 

UTW) 



■FCT> m + m 2 -(l-4/3)(2-/3) 



(A2) 



Hn 



2 

mf , 
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respectively. For the quadratic expression 1A11 in terms of q 2 , we further introduce another determinant Ai such that 



2 , m m (l-FCVn 



Ai = d - 4C2C0 = (1 + 2/3) [3m " 2(1 ~ ~ 4Pmi + 8P)m ~ (1 ~ 4/?)(1 + 4/3 " 8(3 )] • (A3) 

The determinant A of quadratic equation 14811 in terms of D 2 can be either positive or negative in general for different 
parameter combinations of /3, m, q and F. It is fairly straightforward to obtain D 2 solutions by assigning values of all relevant 
parameters. In particular, for one important case of j3 = (i.e., the MSID case with a flat rotation curve and an isothermal 
equation of state), the condition that A > for m > 1 and A = when m — 1 holds for all values of q 2 in a full MSID with 
F — 1 and the condition that A > for m > 1 holds in a partial MSID with < F < 1 (Lou 2002). 

For the magnetic force- free case of j3 = 1/4 with a declining rotation curve with increasing r and a polytropic index 
n = 4/3, we have A > for m > 1 and an arbitrary q > in both full and partial magnetized discs. We now proceed to prove 
this statement rigorously. For /3 = 1/4 of the aligned case, we explicitly have 

c 2 = [m 2 (2FCV m /3 + l)] 2 , ci = 2(1 - FCV m )m 2 [A m {2FCV m /3 + 1) - 4H m /3] , c = A 2 m (l - FCVmf , (A4) 
where < 1 — FCV m < 1 and A m > for f3 = 1/4 and Ai introduced above simply reduces to 

Ai = -24W m (l - FCV m ) 2 m 4 . (A5) 

It is apparent that for Tt m > 0, we have Ai < 0. As C2 > 0, it follows immediately that A > in definition lAlfl . On the 
other hand for 7i m < 0, we have C2 > 0, ci > 0, Co > and therefore A > by definition 1A11 . 

For P = 1/4 in the logarithmic spiral case, we again have the determinant A > for quadratic equation l)69[l . where in 
parallel with definitions 1A11 and 1A21 of the aligned case, the relevant parameters in the spiral case are 



C2 



• C' ) ( ^FCMm + l 



2 , ,2 J . /2„„., . .\ 4 



, ci = 2{l-FCN m ){rn +r) 



A m \ -FCNm + l\--Hr, 



c =A 2 m (l-FCN m ) 2 , (A6) 



with < 1 - FCN m < 1 and A m > for /3 = 1/4. And Ai simply reduces to 

on 

Ai = -yH m (l-FCAT m ) 2 (m 2 +C 2 ) 2 (3C 2 + 27/4) . (A7) 

Hence, as the aligned /3 = 1/4 case, it is apparent that for 7i m > 0, we have Ai < 0. As C2 > 0, it follows immediately that 
A > 0. On the other hand for Ti m < 0, we have C2 > 0, ci > 0, Co > and therefore A > . 



APPENDIX B: PHASE RELATIONSHIPS AMONG PERTURBATION VARIABLES 

In this appendix, we first provide a proof that for the aligned case of (3 = 1/4, the physical D 2 solutions of the plus j/i > and 
the minus 7/2 > determined from quadratic equation 1831 correspond to the factor D 2 — (1 — FCVm) I {FCVm +3/2) being 
positive and negative, respectively. According to Appendix A, the determinant A for both aligned and spiral cases remains 
always positive for the case of /3 — 1/4. 

We start from equation 1831 by casting it in the following compact form of 

P3 — , (HI) 

where y = D 2 and in reference to equation 1831 . pi, p2 and p3 are constants defined explicitly by 

{I- FCVm) {1-FCVm) m 2 q 2 (2FCV m /3 + l) 

S 31 = — I 0^2 olnV77nZ^2 , o/rA ' P 2 = TE^ — TTToT ' 



FCVm + 3(m 2 - 3/2)/(2m 2 + 3/2) ' 5 " (FCV m + 3/2) ' 5a (m 2 + 3/4)[FCV m + 3(m 2 - 3/2)/(2m 2 + 3/2)] ' 

(B2) 

It is apparent that p2 is always positive. While pi and p3 remain always positive for m > 2, their signs depend on the value 
of F for the special case of m = 1. In fact, both pi and p3 diverge at a critical value of F for the m = 1 case as already 
discussed in subsection 3.1.3. 

Let us first analyze the simpler cases of m > 2. As pi > , p2 > , p3 > and pi > p2, we have two D 2 solutions 

pl + P3 + [(pi + pa) 2 - 4p 2 p 3 ] 1/2 ^ P2 + p3 + |p2 - p3 . 

yi = - 2 > 2 > P2 , 

3i + P3 - [(pi + p 3 ) 2 - 4p2p 3 ] 1/2 _ pi + P3 - [(pi + P3 - 2p 2 ) 2 + 4 p 2 (pi - p2)] 1/2 ( B3 ) 



2/2 



2 2 

. Pl + P3 - |pl + P3 - 2p2 
< ~ 2 < P2 • 

We next examine the m = 1 case. The situation becomes somewhat involved since pi < 0, p3 < 0, |pi[ > p2 > when 
< F < F c where the critical F value F c = 3/[7C"Pi(l/4)] ~ 0.6842 is determined from the condition Hi = (see subsection 
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3.1.3). On the other hand, we have pi > p2 > 0, p3 > when F > F c in general. For the aligned case§ and for < F < F c , 
we therefore have two D 2 solutions satisfying inequalities 

Pi + p3 ~ [(pi + P'i f ~ 4p 2 p3] 1/2 „ 

yi = g < o , 

n . , Pl + P3 + [(pi + P'i f ~ 4p 2 p 3 ] 1/2 

< 2/2 = 

2 (B4) 
_ pi + P3 + [(pi + P3 - 2p 2 ) 2 + 4p 2 (pi - p 2 )] 1/2 



. Pl + P3 + |pl + P3 ~ 2p 2 
< ~ 2 = P2 , 

and for F > F c the situation remains the same as <B3t , 

For the spiral case of (3 = 1/4 in parallel, we can repeat this same procedure of analysis. As before, we cast quadratic 
equation I76H in the form of 

(m 2 +Z 2 )q 2 (2FCNm /3 + 1) 

(m 2 + C 2 + 3/4) [FCA^ + 3(m 2 - 3/2)/(2m 2 + 2£ 2 + 3/2)] 

= D 2 {D 2 - (1 - FCN m )/[FCN m + 3(m 2 - 3/2)/(2m 2 + 2f 2 + 3/2)]} 
D 2 — (1 — FCAfm) / (FCAf m + 3/2) 
which yields two distinct D 2 solutions according to the analysis in subsection 3.2.3. The following treatment is identical to 
that for the aligned f3 = 1/4 case described in the first part of this appendix. 



APPENDIX C: ANGULAR MOMENTUM FLUX 



We here provide a proof to show that the angular momentum flux associated with stationary coplanar MHD perturbations 
due to the gravity torque A G , the advective transport A A and the magnetic torque A B remains constant at all radii for the 
entire range — 1/4 < (3 < 1/2. 

We start from expression l|89^ for 5R(</>i ) to derive the following two expressions in a straightforward manner: 

9St(4n) 



89 

8$t(<f> 



2nGm 



p oo p oo 

a / dke~ k ^ z ^J m (kr) / r~ 1/2 sin(m6l - £ In r a )J m {kr a )dr a , 
Jo Jo 



dr 



= -2ttGct 



dk e 



I -k'\z\ 



k'j m+1 (k'r) H J m (k'r) 



Jm(fe'r) 



/ -1/2 



cos(m0 — f \nr' a )Jm(k' r' a )dr' a 



Then in integral iJHUJ, we first integrate the kernel over 9 to yield 



-4tt 3 G 2 - 2 



f oo r oo 

o 2 mr / / F{k,k')e- {k+k ' )M J m {kr)J m {k'r) 
Jo Jo 



k'j m+ i(k'r) H J m (k'r) 



dkdk' , 



where 
F(k,k') 



oo r co 



-1/2 ' -1/2 



J m (kr a )J m (k'r' a ) sin 



On ^ 



dr a dr' a 



We next integrate over the vertical coordinate z to yield the final result 



-47r 3 G 2 <r 2 mr 



Jo Jo JkTkT) Jm{kr)Jm{kr) 



— k'J m +i(k'r) H Jm(k'r) 

r 



dkdk' 



By introducing the following integral transformation 
kr — » :r, k'r — * x , r a —t rrt, r' a —t rr'b 
we reduce expression EH to 



-Air 3 G 2 a 2 m 



where 



OO f oo 



Jm{x^\Jm(x ) 

a: + a; 



x'j m +i(a;') + mJ m (x') 



dxdx' 



M(x,x') = J J r~ 1/2 r' b ~ 1/2 J m (xr b )J m (x'r b )sm fin (^~J 



fin 



(CI) 

(C2) 
(C3) 

(C4) 
(C5) 
(C6) 
(C7) 



where A/" m is the Kalnajs function (Kalnajs 1971). 

With these technical preparations, we finally determine the angular momentum flux associated with coplanar stationary 
MHD perturbations due to the gravity torque as 

(CM) 



AG 2 /~t 2 1 — <G 

— —m-K Ga a 



§ Note these expressions of yi and 1/2 are consistent with the general solution form H51I for the aligned case. 
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Figure CI. Variation of S G versus £, shown in heavy solid line. We plot |0g| 2 in light dashed line for comparison. For increasing £ > 1, 
3p oscillates around zero and its envelope profile fits the |©s| 2 curve with an increasing accuracy. 



F=1, m=2, p=1/4 



\ q=0.3 
\ 



q=1 




F=1, m=2, p=1/4 



\q=0.3 



q-1 



(a) Variation of E. A versus £ 



(b) Variation of S B versus £ 



Figure C2. Variations of S A and H s versus £ for the two branches of y = D 2 solution with m = 2, j3 = 1/4, F = 1, and <j = 0.3, 1. 
Solid and dashed curves stand for y\ and j/2 branches, respectively, with heavy and light curves for q = 0.3 and q = 1, respectively. Note 
that E' 4 is positive and negative for j/2 and y\ branches, respectively, while =. B is positive for both y = D 2 branches. 



where 

I M[x, x) j( x \j( x '\ \— x 'j m+ i (x) + mJm (x'))dxdx' , 
O + x ) 

which is clearly independent of radius r. In the tight- winding regime or the WKBJ limit (i.e., £ S> 1), we may directly apply 
the result of Fan & Lou (1999, see their eqn. 4.13) by noting kr = £ such that 

A G = -sgn(£)m7r 2 Ga 2 |e s | 2 , (C9) 

where 0s has been defined by equation 1811 : here, A G is inward for leading (£ > 0) spiral waves and outward for trailing 
(£ < 0) spiral waves. In general, H G can be calculated via numerical integrations and the specific result is displayed in Fig. 
CI. For fairly open MHD density waves (£ ~ 1) with H G > 0, the angular momentum flux associated with the gravity torque 
is therefore inward for leading (£ > 0) spiral waves and outward for trailing (£ < 0) spiral waves; with increasing f > 1, S G 
falls off rapidly and oscillates around zero with the envelope profile closely fitting the WKBJ result better and better. 

Meanwhile, we compute the other two perturbation contributions to the angular momentum flux, namely, those from the 
advective transfer and the magnetic torque. For the part of advective transport, we start from equation 185H and make use of 
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A A = r 2 S / d0R(ui)R(ji/»O = m ^'" r K|e s | 2 - 9f(8 s )] = m^ 2 Ga 2 H A , (Cll) 



phase relationship 18011 to explicitly derive 

mQ.a?H(e s ) _ 1/2 . mr2cr9(e s ) _ 1/2 

SR(ui) = -r ' sin(£ lnr — mo) H =-= -r cos(£mr — mfl) , 

So So 

M(jiM = -(n<r/E )[l - (l/2 + /3)K(e s ) -C9(e s )]r- 1/2 cos(^lnr -m6») ( C1 °) 

- (n<7/Eo)[(l/2 + /3)9(0 s ) -£»(e s )]r- 1/2 sin(£ In r-mfl) . 

It immediately follows that the advective torque is given by 

r ' 2n , s rmra 2 n 2 r 

where 

= 2 J p 2 (i + 2/3)[c|e s | 2 -a(e s )] 

FC[1 + D 2 - (1 - 4/?)q 2 /(2 + 4/3)] ' 

that is independent of radius r. By specifying parameters j3, m, F, q and £, one finds that for the yi branch £|Osi 2 — 9(Os) < 
for £ > and £|6sj 2 — 9(©s) > for £ < 0. Therefore, the advective angular momentum flux is radially inward for leading 
stationary spiral density waves and radially outward for trailing stationary spiral density waves for the yi branch of D 2 
solutions. For the yi branch in comparison, the situation switches and the advective angular momentum flux is outward for 
leading stationary spiral density waves and inward for trailing stationary spiral density waves. Practically, to achieve the j/2 
branch of stationary logarithmic spiral configurations requires a much lower D 2 which would be unusual for real magnetized 
disc systems. An example of illustration for m = 2, F = 1 and = 1/4 is shown in panel (a) of Fig. C2. 

For the momentum flux associated with the magnetic torque, we start from equation 1871 and make use of phase rela- 
tionship 18UH to first derive 

3?(6 r ) = [mBo<7K(e s )/S ] r ~ 3/2 sin (£ In r - m9) + [mB O 0$S[Qs)/Xo] r~ 3/2 cos(£ In r - m6) , 

tit no (C12) 

3t(bg) = [iB a^{Q s ) /S ] r _3/2 sin(£ In r - m6) + [£B aQ(@ g) /E ] r i/2 cos(£ lnr - m6) . 

It immediately follows that the momentum flux A s associated with the magnetic torque is given by 

a s = r ^ r dzM(b r )U(be) = _ ^yi|e s | 2 = _ m7V 2 Ga 2 E B (C13) 

4tt J J_ x S r 

where 

2g 2 e|Os| 2 



FC[1 + D 2 — (1 — 4/3)q 2 /(2 + 4/3)] ' 

which is again independent of radius r; here, A s is inward for leading spiral density waves (i.e., £ > 0) and outward for trailing 
spiral density waves (i.e., £ < 0). An example of illustration for m = 2, F = 1 and j3 = 1/4 is shown in panel (b) of Fig. C2. 

By comparing the three contributions to the total angular momentum flux transport, we reach the following conclusions. 
First, the contribution from the advective transport part tends to be more dominant over those from the magnetic torque 
and the gravity torque. Secondly, for the y\ branch for stationary fast MHD density waves that has a counterpart in the 
unmagnetized case, leading spiral density waves (£ > 0) transport angular momentum inward and trailing spiral density 
waves (£ < 0) transport angular momentum outward. Thirdly, for the yi branch for stationary slow MHD density waves 
caused by the very presence of the azimuthal magnetic field, leading spiral density waves (£ > 0) carry angular momentum 
outward and trailing spiral density waves (£ < 0) carry angular momentum inward. Note that all these results are only valid 
for stationary logarithmic spiral MHD density waves. There might be some reasons to rule out the yi branch in galactic 
contexts. For example, the yi branch is of lower D 2 which cannot be sustained in real galactic discs. Furthermore, for rotation 
curve parameter /3 7^ 1/4, the yi branch may be ruled out by inequality 1171 . 
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